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We extend the standard intertwining relations used in Supersymmetrical 



^ ' (SUSY) Quantum Mechanics which involve real superpotentials to complex 

superpotentials. This allows to deal with a large class of non-hermitean 
Hamiltonians and to study in general the isospectrality between complex 
potentials. In very specific cases we can construct in a natural way "quasi- 



■ complex" potentials which we define as complex potentials having a global 

^ : 

, property such as to lead to a Hamiltonian with real spectrum. We also 

OO 

0^ I obtained a class of complex transparent potentials whose Hamiltonian can be 

^ • intertwined to a free Hamiltonian. We provide a variety of examples both for 

. the radial problem (half axis) and for the standard one-dimensional problem 

a ; 

^ ' (the whole axis), including remarks concerning scattering problems. 



1. Introduction and general formalism 

Non-hermitean interactions have been discussed in Field Theory and Statistical 
Mechanics, see for instance with applications to Condensed Matter and to 
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Biology 0. In the context of Hadronic and Nuclear Physics, the so-called optical 
potential framework also involves this type of Hamiltonians as studied by Baye et 
al. P]: the imaginary part in the potential signals the opening of new channels. In 
the present paper we investigate this class of systems within the Supersymmetrical 
(SUSY) approach in Quantum Mechanics with particular attention to complex po- 
tentials associated to a real spectrum. Recently there was a lively activity on complex 
potentials with the property of PT (parity and time reversal) invariance, associated 
to real and positive spectra 

In the standard approach to SUSY Quantum Mechanics one makes use of inter- 
twining relations where the intertwined Hamiltonians are hermitean and the inter- 
twining operators depend on real " superpotentials" 0. 

Here we generalize this approach to complex superpotentials W{y) which allow 
to deal with non-hermitean Hamiltonians. Our approach makes use of supercharges 
with complex superpotentials intertwining two partner Hamiltonians [§: 

= (1) 

The Hamiltonians H^^\H'^^'' contain complex potentials V^^\y) = U^\y) + iUf\y) 
and V^^\y) = U^\y) + iUj^\y) and the supercharge is defined by 

q+ = -d + W{y) = -d + f{y) + ig{y) (2) 

with f{y),g{y) real functions. The role of the intertwining operator g"*" is to map so- 
lutions of the Schroedinger equation corresponding to H^^^ to solutions corresponding 
to H^^^ for identical spectral parameters \1/ '•''-'(?/) ~ q^'^^^^y). 

By comparing (|l]) with the factorization approach 0, [^] one can explicitly derive 
that (0) is equivalent to the factorization at a complex energy E = + iej of the 

above mentioned Hamiltonians 0. Furthermore by taking successively the hermitean 
"'We know however that intertwining is more general than factorization as one can check from 
matrix SUSY Quantum Mechanics |^ and from Higher order SUSY Quantum Mechanics 
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and complex conjugations a second intertwining relation can be derived: 

g-iy(o) = H^^)q-. (3) 

The supercharge 

q- =d + W{y) = d + f{y) + ig{y) (4) 

is related to by successive hermitean and complex conjugations. It maps solutions 
of the Schroedinger equation corresponding to H^^^ to solutions corresponding to H^^^ 
for identical spectral parameters '^^^\y) ~ 

In general one expects that these operators will not map physical wave functions 
into physical ones. We will study the conditions allowing to realize such a mapping 
apart from one state, solution of the equation = Q , called zero mode : from now 
on we will refer to such a mapping as to the case of strict isospectrality. 

In the standard factorization approach to SUSY Quantum Mechanics was suit- 
ably extended to the case of complex potentials. In the construction of isospectral 
Hamiltonians our emphasis will be more on the role of the intertwining relations (0) 
with complex superpotentials W{y) = f{y) +ig{y), including the case of intertwining 
hermitean and non-hermitean Hamiltonians. 

Eqs.(|l]) and (|^) yield the following relations between the potentials V^'^\ V^^"^ and 



the super potential W{y) : 

uPiy) = -ny) + f{y)-g\y) + en; (5) 

U?\y) = -9'{y) + 2f{y)g{y) + ej; (6) 

u^R\y) = f\y) + ny)-9\y) + en, (7) 

u^Piy) = 9\y) + V{y)9{y) + ei, (8) 



where enj are real numbers which in the factorization approach correspond to the 
factorization energy E = en + iej. 
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These relations express both complex potentials in terms of two arbitrary real 
functions f{y),g{y) determining the complex superpotential 0. As far as their solution 
is concerned we would like to stress that in general it is not possible to consider one 
potential to be given and the other to be constructed since this problem in general 
cannot be solved analytically. What one can do is to provide interesting classes of 
model systems with Uj^\ Uj^\ f,g to be determined contextually. 

The problem can be formulated both as one dimensional (the whole axis, y = x) 
for the wave function \E^(x) or as a radial problem (half-axis, y = r) for r-ip{r) = \['(r). 
The conditions for mapping physical wave functions into physical wave functions can 
be different for these two cases. 

There is a possibility of imposing U^\y) = uf'\y) or U^^\y) = U^\y) or more 
complicated conditions since we deal with two functions /, g. The first conditions 
leads to g' = but not necessarily to trivial Uj whereas the second condition leads 
to /' = 0. For the whole axis problem both conditions allow to map physical wave 
functions to physical wave functions. For the radial problem the bound state wave 
function r ■ ilj{r) should vanish at the origin like r'"*"^ and decrease at infinity in order 
to be regular and normalizable. Regularity and normalizability should be preserved 
in the mapping between physical wave functions: while g{r) can be assumed to be 
regular at the origin, this implies for the function /(r) at the origin: 

/(r) or f{r) ~ ^. 

In the factorization approach W{r) = — ^^■''/\&^'* , where \E'^'* is a solution (physical 
or unphysical) of the Schroedinger equation for E = + iej : its behaviour r'^^ 
(physical) or r~' (unphysical) corresponds to the two behaviours of /(r), respectively. 
This expression for the superpotential W{r) displays its parametric dependence on 
E = en + iej. The first behaviour maps the centrifugal barrier for I to the partner 
^Indeed V^^\y) - V^°^{y) = [q-,q+] = 2W'{y). 
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barrier / + 1 while the second leads to / — 1. So for the radial case the condition 
/ = const is not compatible with the physical behaviour at the origin. 

In contradistinction, on the line the boundary conditions for the wave functions 
are a less severe constraint because one can have arbitrary smooth behaviour at the 
origin and suitable decrease at ±oo. 

In Sect. 2 we derive, from the basic formalism for the construction of a SUSY 
partnership for complex potentials given in (|^) - (|^), the conditions which allow 
for a partnership between a hermitean Hamiltonian and a non-hermitean one. The 
discussion is formulated for the half line (radial problem). 

Section 3 concerns with the full line problem (one-dimensional Quantum Mechan- 
ics). We construct complex transparent potentials (Subsect.3.1) and give examples 
of non-polynomial "quasi-complex" Hamiltonians with real spectrum (Subsect.3.2) 
including remarks about their symmetry properties . 

An outlook on the general isospectrality between complex potentials is given in 
Sect. 4 with contact to the Zakharov-Shabat problem, to scattering problems and to 
Higher order derivative SUSY transformations. 

2. Uj^\r) = and "quasi- complex" potentials for the radial 
problem. 

We study a special case of partner relationship in which one partner is real and 

the other is complex. The non obvious aspect of this relation lies in the fact that 

the spectrum of the complex potential will be real only if we are able to satisfy the 
^As mentioned in |^ there is an alternative interpretation which considers stiU the angular mo- 
mentum to remain I but allows for the existence of a repulsive centrifugal like singular potential of 
dynamical origin. 
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physical conditions at the origin discussed above. Otherwise the operator intertwining 
relations (|l]), @) involve unphysical solutions of the Schroedinger equation. 

We try to find cases where the real potential U^\r) has many bound states be- 
cause it is not astonishing that in the scattering regime there is a real continuum 
spectrum for complex potentials as one can verify in the case of a constant complex 
potential (cf. also eq.(39) of and its subsequent complex extension). Note inci- 
dentally that, at least restricting to real potentials, the constant potential and the 
example of |^ are well known standard SUSY partners. 

Eq.(|^) for Uj'^\r) = allows to express the function /(r) in terms of g{r) : 



Thus we can obtain a complex potential V^^^ which is equivalent, in the sense of 
intertwining relations, to a purely real potential. From now on we will call such 
complex potentials "quasi-complex" potentials. The formal definition is given in 
terms of the equations above. 

The requirement that fiux is absorbed and not generated leads to the dissipative 
condition Uj^\r) < i.e. g'{r) < 0, which is quite common in the context of optical 
potential studies. In this case the previous equations cannot be solved unless U^^^ (r) 
has no bound states. This is so because the eigenvalues of the Hamiltonian for a 




(9) 



which implies: 




(10) 




(12) 



7 



complex dissipative potential have negative imaginary parts [] 

, P [2!,-g'{r)\%Xr)?dr - f (^Ur)^n(r) - ^n(r)vl/;(r))-] 

(Jo l^n(r)pdr) 

while the partner (hermitean) Hamiltonian should have the real spectrum. 

It is instructive to check how this is realized for a class of models for which the real 
part /(r) of the superpotential coincides with the class of solvable models in standard 
SUSY Quantum Mechanics [PH]]: 

/(r) = A ■ tank ar — B ■ coth ar. 

One finds that the condition of dissipativity {g'{r) < 0) implies g'{0) ■ g'ioo) > and 
leads to the requirement S < 0, in conflict with the condition A > B > required 
for the existence of bound states for the potential U^^\r). 

Not imposing the dissipative condition we provide a variety of examples with a 
pure discrete spectra. The centrifugal behaviour of /(r) at the origin implies (§) for 

Jv)^ ^ gix) a ■ T] a 



2/o + 3' 



/(r) ~ - g{r)r^a-r; 



r ' ' 2/o - 1 

or for 6/ = : 

r 

Given this behaviour at the origin for e/ 7^ 0, it is natural to solve in terms of the 
ansatz: 

g{r) = a ■ r(j){r); 0(0) = 1; 0(oo) = 0, 

while the first condition for 0(r) is obvious, the second one is suggested by the re- 
quirement of having a purely discrete spectrum (|lO]),(|lT]). 

^The terms not appearing in eq.(8) of jsj vanish if one imposes physical boundary conditions for 
the wave functions. 
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First we construct potentials which have the behaviour of generahzed oscillators at 
infinity with the mapping which preserves normalizability and therefore corresponds 
to strict isospectrality. We can choose among a variety of cases 

0(r) = (r'- + l)-i; m > 1, (14) 

Asymptotically one obtains for large r: 

/(r)^^-r^"^-^ U^R\r)-f; ^'^^\t) ^ expi-pr'^)- (15) 

with p a real parameter. 

Alternatively we can choose: 

0(r) = exp{—ar)] o; > 0. 

Asymptotically for large values of r: 

/(r) ^ ^ • expiar)- t/f (r) ~ /^(r); (r) ~ exp{-p ^^^dr); 

Despite the exponential growth of /(r) one can check that the mapping preserves 
normalizability because of the extremely fast decrease of the wave functions. 

3. uf\x) = : transparent and "quasi-complex" potentials 
on the line. 



In Sect. 2 we have discussed intertwining relations between real and complex po- 
tentials in SUSY Quantum Mechanics which allowed us to introduce for the half line, 
i.e. the radial problem, the concept of quasi-complex potentials. Due to the less strin- 
gent conditions at the origin we will be able on the line to introduce within the class 
of quasi-complex potentials a sub-class of complex transparent potentials, equivalent 
in the sense of the intertwining relations to a vanishing potential. 
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3.1. Complex transparent potentials on the line. 

We begin from a specific attempt, the so called transparency problem (cf. for 
example 2.5 of to intertwine a complex Hamiltonian H^^'^ to a free Hamiltonian 
with f/}j°^ = uf^ = 0. We insert the condition e^j = e/ = in eq. (|TUp and 
impose its vanishing: 

Making the substitution g{x) = ±p^'^{x) one obtains 



which can be integrated and yields: 

a —1 

9{X) = 57 TTTT = Im ; — (17) 

^ ^ 1 + a^{x + b)^ x + b + i/a ^ ^ 

in terms of the arbitrary real constants a ^ and b. By inserting (|T^ in ( ]TT| ) we are 
thus able to write the "transparent" complex potential 

V^^\x) = 2a^tl±^^t±3_,, 4a3(. + 6) 2 



[1 + a2(x + 6)2]2 [1 + a2(a. + 5)2]2 (3. + 5 + 1)2- v j 
This potential gives a trivial S-matrix but generates a zero-energy bound state, 

*o'^(^) = — ?-T- (19) 
x + b+ ^ 

The potential V^^\x) is invariant under PT reflection where the P-parity inversion 
is performed with respect to x = — 6 and the T-transformation leads to the complex 
conjugation of the potential ( see discussion on the role of PT-invariance in next 
Subsection, after (p7|) ), the imaginary part of the potential is odd under P-inversion. 

The transparent potential (0) is a limiting case {en 0^) of the analytic solution 
of ( p!OD with e/j < and still ej = 0. Proceeding as before we can integrate 

{p'\x) + p-\x) + eRp\x)y = 
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and obtain (with a, b real constants) 



g{x) = ±p-^{x) = ± — —— -, (20) 

a — va^ — Aeji cosh[2^/ —eji{x + b)] 

where g{x) has in fact no singularity for real x because is negative. In the limit 
0^ one can indeed check that ( pOD reduces to (p!7[) . Without giving the details 
of the derivation we write the resulting complex transparent potential: 



(a^ - ien) - a J (a^ - AtR) cosh{2^-eR{x + b)) 
U\,\x)=AeR- 



^{a-' - AtR) cosh{2^^{x + b))Y 



Ui^\x) = ,^^V^n^(-'-^^n)s^nhi2V^ix + b))_ ^^^^ 
[a - ^{a? - AeR) cosh{2^^{x + b))]^ 

The total superpotential W{x) can be reexpressed as: 



W{x) = -V^-tanh{y/^{x + b) +ip), (22) 

if instead of a we introduce a real constant p = ±1/2 ■ Arctan{2y/—eR/a) with p ^ 
|(2n + 1) in order to avoid potentials with strong singularities for real values of x. 
Correspondingly, (0) can be simplified: 

V^'\x) = U^^\x)+tUi'\x) = — . (23) 

^ ^ ^ ^ ^ ^ ^ ^ cosh^{^/^{x + b)+ip) ^ ' 

One can easily check that it vanishes at infinity and has a bound state with energy 
eR given by 

^^^i(^)= , V (24) 

cosh[y'-eR[x + bj + ip) 

The potential is invariant under PT transformation and its imaginary part is P-odd. 

It is worth to note that there is here a very specific reason why the complex poten- 
tial can be obtained from well-known real transparent potential cosh^"^ by a complex 
shift of the coordinate x. Indeed, the Hamiltonian H^'^^ with constant potential is in- 
variant under an arbitrary (even complex) translation in x. Factorizing H^^\ we can 
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thus obtain a variety of superpotentials which differ by this arbitrary shift of x. Cor- 
respondingly, we have a variety of superpartner Hamiltonians H^^^ with potentials of 
the form (^31) both with the real and complex shifts of x. In the framework of SUSY 
transformations, Hamiltonians of this variety can be connected each to other by the 
second order SUSY transformations via the intermediate Hamiltonian H^^\ 

It is necessary to stress that in general a complex shift of x can lead to drastic 
and non trivial consequences. In particular, it can change the analytical properties 
of the potential (and of the wave functions): compare (for real x) the nonsingular 
potential (|l^) with its very singular real analogue 2/(x + 6)^. On the contrary, SUSY 
transformations with nonsingular superpotentials represent the tool to build com- 
plex potentials with spectral properties under control (mapping eigenf unctions into 
eigenfunctions). 

Transparent potential for e^j > and e/ = can be formally deduced from ( pOD 
with hyperbolic cosine appropriately replaced by trigonometric cosine. Again there 
are no singularities in g{x). The resulting potential has oscillatory behaviour, ap- 
parently compatible with transparency, taking the limiting case eR 0^ leads to 

3.2. Non-polynomial quasi-complex potentials on the line. 

In this Subsection we construct quasi-complex potentials on the line restricting 
ourselves to potentials bounded from below at infinity f\. Let us first illustrate why we 
are unable to construct quasi-complex polynomial Hamiltonians by the intertwining 
relations (|I]). If we assume that g{x) ~ x" and we require the potentials to be 
bounded from below, from Eqs.(^ - (|^) we can derive that f{x) grows at oo as a 
^Scc [T^ for discussion of eigenvalue problem for real potentials unbounded below at infinity. 
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power or faster. This excludes the possibihty to obtain uf\x) = with this ansatz, 
at least with the supercharges (|^), since eq.(^ is not compatible with the growth 
of f{x) at oo. For completeness we mention that the cases of constant f{x) or g{x) 
are presented in Sect. 4 and that some remarks concerning non-hermitean polynomial 
Hamiltonians will be given at the end of this Subsection. 

As a first class of examples of non-polynomial quasi-complex potentials it is possi- 
ble to extend the generalized oscillator ([1^), (|15D to the whole axis. As for the radial 
case one obtains strict isospectrality since the mapping preserves normalizability. 



For real potentials it is well known [10| that solvable models can be generated by 



assuming a superpotential W{x) ~ tank ax. Here we investigate the ansatz for g{x) 
to be inserted in (|rUI),(|TTl),(|g): 

g{x) = P ■ tanh{ax). (25) 



This leads naturally to a definite sign for g'jx), cf. (pr2D, for example dissipativity if 
aP < 0. iFrom eq.(|D one can derive that /(x) at ±oo becomes constant. 

Imposing furthermore the absence of centrifugal like singularities at the origin, 
implies a numerical condition on e/, namely ei = a ■ f3. 

From general arguments dealing with a dissipative case (|T3|) strict isospectrality 
implies that there cannot exist bound states if the mapping preserves normalizability. 
Indeed from isospectrality, one would deduce that in the same time eigenvalues have 
imaginary part for one potential and are real for the other, but they should also be 
equal. Of course all our arguments are to be intended to be true in general apart 
from the possible existence of zero modes. 

First we argue that no bound state exists for E = + iej. This can be derived 
by taking into account that the zero modes of the operator q" (^) correspond to the 
solutions of the Schroedinger equation for complex spectral parameter E = + iei. 
From this one can conclude that the asymptotic behaviour of f{x) is not compatible 
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with the existence of a normahzable solution for that energy. 

The absence of bound states at other energies is however now not necessarily 
paradoxical since the potential is finite. 

Although the previous arguments provide a proof that there cannot be bound 
states it is useful to develop an independent proof of the absence of bound states for 
the real potential 



4^2 2g V 



(^) = 7^ - ^ - ^ + - + ^R- (26) 



The task of a proof of absence of bound states is normally not an easy one for a 
generic potential [13|. In our case, however, performing the substitution of (|25|) in 
(p6|) one can directly conclude that the potential is strictly repulsive (independently 
of the value of a and f3) and proportional to sech'^ax, apart from an energy shift. A 
discussion of this potential in the radial case can be found in 0. In contradistiction 
U^\x) can be attractive. In particular one can choose the constants a = —2(3/y/3 
such as to make Uj^\x) =const, which renormalizes the real part of the energy. This 
last choice leads to an interesting relation between a genuine repulsive real potential 
and an absorptive purely imaginary potential. 

In this Section and in Sect. 2 we have constructed quasi-complex potentials which 
have nontrivial imaginary part but the spectra of the corresponding Hamiltonians are 
real apart from (possibly) one level. From eq.(|I^) it is clear that in some peculiar 
average sense the imaginary part Uj^\x) of the potential vanishes: we stress that 
Eq.(p!3D vanishes for all eigenfunctions, thereby reflecting a global property of the 
potential 0. 

We give a more technical illustration of such a type of global property and discuss 
^It is useful to remark that introducing explicitly h in the Schroedinger equation and taking the 
classical limit h the Eq.(|l2|) vanishes and (|lO|) coincides with (jl^) because of the last three 
terms. 
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how it can be implemented on the full line by a suitable Gaussian ansatz with e/ = 
: g{x) = 7 ■ exp(— ax^") with positive a and 7, qualitatively to be interpreted in 
terms of "range" and "strength" parameters. Correspondingly f{x) = g'{x)/2g{x) = 
— arax^""^. By this ansatz we construct the real potential V^^\x) and the complex 
potential V^^\x) according to ([10|) - ([121) . 



The potentials one obtains are: 

(x) = f/l*) (x) = a^n^x^"-^ + an{2n - l)^'""' - 7' ■ exp(-2«x'") + e^; 
(x) = (x) + 2f/f ^ (x) = aWx^""-^ - an(2n - l)a;='"-2 _ 

- 7^ ■ exp(-2ax^") + e/j - 4m7ra^'""^exp(-ax^"). (27) 

Both potentials are regular and growing at infinity leading to a discrete real spectrum 
since the supercharges respect in this case the normalizability of the wave functions. 

We make explicitly use of the supersymmetry to connect, apart from a normal- 
ization factor, the eigenf unctions '^'^^^{x) and \E''^^)(a;) and rewrite (|T^) 

ImE^ = ^ ^ " ^ ^' 28 

with 

¥^\x) ~ (9 + fix) + z^(x))^(°)(x). (29) 

With our choices IniEn = because |\I'^^-'(x)p is even under parity while g'{x) is 
odd. Indeed H^^^ is invariant separately for parity P and time reversal (complex 
conjugation) T inversions while H^^^ is invariant only under the combined PT trans- 
formation. So ^'^'^''(x) has definite parity and can be chosen real whereas '^^^\x) is 
complex and has no definite parity, but by the previous expression (|29|) we control the 
parity invariance of |\l/^^'*(x)p. Similar arguments hold also for the case of transparent 
complex potentials ([T^), (PB|). 

Heuristically one can speculate that in some suitable limit of a small and 7 large 
the purely polynomial terms can be neglected. Restricting to a finite interval of x, the 
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exponentials can be, roughly speaking, approximated by unity. Apart from an energy 
downward shift, taking such a limit for 0,7 in a finite interval of x, the dominating 
term is uj^\x) = — 4ia7nx^"~^. 

The invariance properties of the model (0) are shared by systems with polynomial 
potentials V = Ax"^^ + iBx^"'^^ on the whole axis. The Hamiltonian H is as before 
not parity invariant and not time reversal (complex conjugation) invariant but still 
possesses an antiunitary PT invariance. This property implies the reality of the 
matrix elements of H^^^ in a suitable basis (for a constructive proof cf. IQ). This 
obviously is consistent with the vanishing of (|28|). 

It is remarkable that this type of systems was discussed recently in a framework 
of a supersymmetric quantum field theory [|T3| (in connection with the studies of the 
so-called Lee- Yang zeros) and the spectrum was found numerically (Blencowe et al. 
IIT5I1) to be real and bounded below. 



4. Outlook and perspectives. 



In the previous Sections we have illustrated the consequences of the intertwining 
relations by discussing several examples in considerable detail. In this last Section we 
want to mention briefly some topics which can also be treated in the formalism we 
have discussed above. For conciseness we will not illustrate in detail various possible 
applications but only outline few paths which can be explored. 

4.1. Scattering. 



Considering potentials with continuum spectrum (see for example [0), we start 
from asymptotic states of a real potential uf\r) = 0. These scattering states of the 
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real potential {r) in partial wave lo, asymptotically at oo (for the notations cf. 
§) read: 

^(°)(r) ~ exp{-ikr) + {-iy''+^S}^\k)exp{ikr). (30) 

We illustrate for a specific value of = fco, how one can generate for the complex 
partner Hamiltonian a partner wave function '^^^\r) with no ingoing wave , a condi- 
tion commonly related to bound or quasi-stationary states (resonances) in the case of 
real potentials. We choose g{oo) = + a ■ exp{—[3r). Assuming e/ = 0, from (P) we 
obtain /(oo) = — (a;/3/2fco) ■ exp{—(3r). By applying g~(oo) = {d + f{oo) + ig{oo)) to 
the ingoing wave in \E'''^^(r) is exponentially damped. The energy associated 
to dH) for A; = A;o is indeed E{ko) = k^ + U^^'^ (oo) = en, as expected from standard 
SUSY Quantum Mechanics, cf. Andrianov et al. in [0. A similar treatment for 
k = —ko leads to an exponential damping of the outgoing wave in '^^^\r) which can 
be interpreted as black sphere or strong absorption limit for this specific value of k. 

Even if in the present paper the main attention has been given to the case of 
strict isospectrality which respects normalizability, there are other interesting cases 
where intertwining relations do not respect normalizability because of the behaviour 
of / and g. For a suitable falloff of the complex potential, bound and quasi-stationary 
wave functions are for large values of r, correspondingly, the decreasing and increas- 
ing solutions of the (asymptotic) equation — ^" = {Eji + lEi)^. One could consider, 
for example, the SUSY partnership between a complex potential, leading to a bound 
state with positive real energy and negative imaginary part of energy and an- 
other complex potential for which at this complex energy one has a non-normalizable 
quasi-stationary (resonant) wave function [|l^. This can be achieved with the same 
techniques outlined above which produce an exponentially decreasing state by ap- 
plication of a suitable supercharge g to a non-normalizable state (creation of bound 
state). Similar results can be obtained by complex rotation approach (for example. 
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T8| and references quoted therein). Alternatively, the reverse mapping from normal- 
izable to non-normalizable (cancellation of bound state) can be realized because (§) 
does not hold for the case of complexity of both potentials. 

For scattering on the line one can construct SUSY transformations for e/ = which 
also lead to an extinction of the ingoing wave. On the other hand transformations 
with 6/ 7^ (an explicit realization can be obtained for the potential defined by 
formally provide the extinction only for complex values of k. 

4.2. Isospectrality between complex potentials. 



This topic has been thoroughly investigated by Baye et al. with special emphasis 
on the construction of phase equivalent potentials. In our approach starting from the 
equations (|^) - (|) one can construct families of isospectral complex potentials if f{y) 
and g{y) are such as to preserve the normalizability in the mapping (cf. discussion 
in Sect. 2 and Sect. 3 of this paper and the approach of (P). This allows in principle 
to investigate the interplay between repulsion /^(y) (cf. eq.(pTD below) in the real 
part of the potential and the absorptive potential (depending on f{y) both for bound 
states and the scattering problem). 

For simplicity we consider the case of the whole axis which e.g. is appropriate 
for physical applications concerning solitons (cf. the so-called Zakharov - Shabat 
problem and its generalizations in 3.9 of [|ll|])- The most simple example has already 
been mentioned in Sect.l and can be formulated as follows: 

g[x) = K, = real const; 
uP{x) = ui'\x) = 2Kf{x) + ej; 
4°)'(^)(a;) = Tf'ix) + fix) + en- k\ (31) 

As a particular case we can solve the above equation for U^^\x) = or U^^\x) = 0, 
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this is the well known analytically solvable case of the Riccati equation ||19|| . 

We can also discuss a mapping|3 which does not affect the real part of the potential 
with f{x) = A = real const: 

U^^\x) = u'i\x) = X'-g'{x)+en; 

uf'>''^'\x) = T9'{x) + 2Xgix) + ei. (32) 

One can notice that for any g{x) , apart from an energy shift, U^\x) = U^\x) are 
attractive. 

Ansaetze leading to rather simple expressions for eq.(^) are g{x) ~ sechx or 
g{x) ~ [2 cosh X ■ seek 2x — 1]. 

Applications involving A = given in include problems with two-state systems 
and some phenomenological models for quark-gluon plasma physics and neutrino 
propagation in matter. One can also notice that the Zakharov-Shabat problem, cf. 



2.12 of 11 : 



n[{x) + iu{x)ni{x) = (n2{x); n'2{x) — iu{x)n2{x) = —(ni{x) 

reduces to the Schroedinger equation with V^^\x) of eq.(|32D for A = : 

—n'l — {u^ + iu')ni = C^^^i- 

In this respect we would like to remark that in the mathematical literature there are 
interesting problems |ll| (Sect. 1.5 and Sect. 5.1) where the imaginary part Ui{x) of 



complex potential is not everywhere dissipative. 
As a typical example. 



Ui(x) ~ (l + 2tanh2x) 

cosh 2x 



■"in Sect. 2 we have explained why this cannot be done on the half axis. 
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results from a function g{x) ~ exp{—x) ■ sech2x with /(x) = 0. In this case Uj (x) 
and Uj^\x) have opposite sign apart possibly for a constant ej and instead Ujl\x) 
and U^\x) are equal. 

As another exercise it is possible to construct the potentials if f{x) in ([3T| ) and 
g{x) in ( |32D are assumed to be of Gaussian shape, like in Subsect. 3.2. 

4.3. Higher order derivative supercharges. 

By suitably iterating intertwining relations with different / and g which strictly 
preserve isospectrality, as explained above, one can generate an infinite class of quasi- 
complex potentials. Finally, one can also find an intertwining relations where the 
Hamiltonian is again hermitean: this provides an explicit construction of a 

Higher order irreducible SUSY transformations as discussed in 0], i.e. transforma- 
tions which cannot be represented as a sequence of two standard (with real superpo- 
tentials) first order transformations. The concept of irreducibility depends on the class 
of intermediate Hamiltonians which are allowed. What was an irreducible transfor- 
mation allowing only hermitean Hamiltonians P] can become a transformation which 
effectively is reducible allowing non-hermitean intermediate Hamiltonians. In a sim- 
ilar way one can discuss second order complex SUSY transformations between two 
non-hermitean dissipative Hamiltonians. Again the intermediate Hamiltonian is not 
necessarily dissipative everywhere. 
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